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Abstract 

Order parameter equations, such as the complex Swift-Hohenberg (CSH) equation, offer a sim- 
plified and universal description that hold close to an instability threshold. The universality of the 
description refers to the fact that the same kind of instability produces the same order parameter 
equation. In the case of lasers, the instability usually corresponds to the emitting threshold, and 
the CSH equation can be obtained from the Maxwell-Bloch (MB) equations for a class C laser 
with small detuning. In this paper we numerically check the validity of the CSH equation as an 
approximation of the MB equations, taking into account that its terms are of different asymptotic 
order, and that, despite of having been systematically overlooked in the literature, this fact is 
essential in order to correctly capture the weakly nonlinear dynamics of the MB. The approximate 
distance to threshold range for which the CSH equation holds is also estimated. 

PACS numbers: 42.65. Sf,05.45.-a 
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I. INTRODUCTION 




The complex Swift-Hohenberg (CSH) equation is an order parameter equation that pro- 
vides a reduced description of a variety of systems such as Rayleigh-Benard convection 
{2! , optical parametric oscillators 0, I^, [^l , Couette flow [gI , nematic liquid cr j 
netoconvection [8] , propagating flame front [9] and photorefractive oscillator 
others. In the case of lasers operating near peak gain (small detuning), a derivation of the 
CSH equation for class A and C lasers was obtained in [12I, 13], starting from the semiclas- 
sical Maxwell Bloch (MB) equations Q, Q, 0, 0, Q, that provide a general description 
of transverse patterns in two levels, wide aperture and single longitudinal mode lasers. For 
class B lasers, such as CO2 and semiconductor laser, CSH equations have been obtained in 
H and y respectively. (For an explanation of the classiflcation of lasers, see Ref. j2l|.) 
2g)erimental observations of patterns in wide aperture lasers were reported in, for example. 
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261]. (For reviews on pattern formation in nonlinear optical systems, see 



) 

The deduction of a generic order parameter equation greatly simplifies the theoretical 
description of the system. But it is important to note the limitation of this kind of model 
equations, as Cross and Hohenberg state in their review [1, p. 874]: "It is true that many 
properties of nonequilibrium systems are encountered in these equations, and indeed many 
hard problems (...) may profitably be addressed in the simple framework provided by these 
equations. However, it is only as a perturbative expansion valid in a small region near 
threshold that they provide a quantitative description of real experimental systems, and 
results may be even qualitatively misleading if applied far from threshold." 

We will focus our attention in a class C laser near peak gain that is assumed to be well 
described by the corresponding Maxwell Bloch equations. Using the assumptions that the 
amplitudes of the physical fields are small and depend slowly on time and on the transversal 
spatial scales, a CSH equation is derived form the Maxwell Bloch equations. The resulting 

n n 

CSH equation is that derived by Lega et. al. in |12L I13I ]. Although not explic- 
itly stated in this original formulation, the equation includes terms of different 
asymptotic order, in contrast to the Ginzburg Landau equation that is obtained 
for negative detuning jsoj]. This asymptotic nonuniformity is systematically obviated 
in the literature and was only recently addressed [si)]. It is the manifestation of the fact that 
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dispersion and diffusion have necessarily different asymptotic order, and it affects 
the slow scales that the system develops near threshold. 

The argument of a quahtative only scope of the model equation is usually invoked to 
justify the application of the CSH equation far from threshold. The qualitative correctness 
is difficult to be theoretically established, but, on the other hand, the capacity to produce 
quantitative predictions can be determined from the numerical integration of both, the origi- 
nal system of Maxwell Bloch equations and the CSH equation. The results of the comparison 
between the numerical simulations of the CSH and the MB equations is what we present in 
the subsequent sections of this paper. We can obtain from the simulations the relative error 
that introduces the approximation and estimate how far from threshold we can increase the 
pump while keeping a small relative error. Also, this numerical comparison between the 



CSH and the MB equations provides a confirmation of the main result presented in 



that the CSH equation for the description of the weakly nonlinear dynamics 
of the system near threshold (derived in 1 I12I . fisl ]) necessarily contains terms of 
different asymptotic order. 



II. THE COMPLEX SWIFT-HOHENBERG EQUATION AND ITS NUMERICAL 
PHASE DIAGRAM 

The Maxwell-Bloch equations for a two-level single longitudinal mode laser with flat 
mirrors are 

^ = iaV^E -aE + aP, (1) 

§ = -{l + in)P+{r-N)E, (2) 

^ = -bN + ^(EP + EP), (3) 

where E{x^ t) and t) represent the complex electric and polarization flelds, and 

N{x^y^t) is the real valued fleld of the population inversion (the same nondimensional 
formulation as in Ref . js^ is used) . Parameter a > is the strength of the diffraction (that 
we set to 1 by scaling the space variables), a > is the cavity losses, is the cavity detuning 
(the difference between atomic and resonance frequencies), r is the pumping parameter, 6 > 
is the decay rate of the population inversion, = / dx^ + / dy^ is the laplacian operator 
in the plane transverse to light propagation, and the bar stands for the complex conjugate. 
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We will consider, specific CB^se, ^ class C laser, for which a ^ 1 and 6^1. 

The corresponding CSH equation was obtained in [I2I, and a simpler derivation 
method, in which no a priori relative scaling of the variables is assumed, was introduced 



m |3l|]. A linear stability analysis of the Maxwell Bloch equations shows that the lasing 
instability takes place at a critical value of the pump Tc = 1. The assumptions of small 
detuning and small distance to threshold, 

<l,|r-l| <1, 

are expressed through a small parameter < ^ <C 1: 

il —CJS, 

a 

r-l = ^^^(cu2 + a)£^ (4) 



were a and u are order 1 parameters that represent the scaled pump and detuning respec- 
tively. 

The resulting CSH equation is of the from: 

= a(^ + iV2(/) - (t)\(t)\^ - 2^cjVV - ^^VV, (5) 

where time and space were scaled as t = ^^^^te'^ and (x, y) = (x, y)£. 

This CSH equation is exactly the same as that obtained by Lega et. al. in 
{12! , 13], but with the variables rescaled to show that it has terms of different 
asymptotic order and that it is not possible to remove the small parameter s 
from the equation. This asymptotic nonuniformity comes from the simple fact 
that dispersion involves second order spatial derivatives while double diffusion 
has fourth order ones and thus, in the long wave approximation where higher 
derivatives correspond to smaller terms, these two terms have necessarily dif- 
ferent asymptotic order. This crucial fact is precisely what forced Lega et. al. 
I2I , [13 ] to derive the CSH expanding first up to two orders (the first one in- 
cluded dispersion and the next the double-diffusion) and then collapsing back 
the expansion to get the CSH equation. But, despite of the wide use of the CSH 
equation, the asymptotic nonuniformity is never mentioned in the literature. 



and it was only recently analyzed in 



3 if ] where it was shown that it gives rise to 



two characteristic slow scales: one associated with dispersion 5disp ^nd a second 



one associated with diffusion 5diff- Using the scahng indicated above, 5disp ^ 1 
and Sam ^ <^ 1? but in the original scahng of the MaxweU Bloch equation 
^^disp ^ 1/s ^ 1 and Sdm ^ Vv^ 1, so both are long spatial scales. 

The CSH equation above has to be considered in the close-to-threshold limit of 6: ^ 0, 
and the relation between the Maxwell Bloch and CSH solutions can be written as 



E{x,y,t) 




1 


P{x,y,t) 




1 


N{x,y,t) 








Vb 



(^ + 1), 



-iawi/e 



a 



ix,y,i)e + 0{s^), 



(6) 



Traveling wave solutions of the form (/)tw = ^/ae^^^'^'^ i^tw^^ with /ctw = Ikxwl ^ 1 are 
approximate solutions of the CSH equation up to 0(e) corrections (this family of TW is just 



the result of making; the 



TW family, see 



12. 
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imit £ ^ and k ^ 1 in the well known expression of the exact 



3, 



33[|). This solution exists only for a > and a linear stability 
analysis shows that it becomes unstable outside the region defined by > and a > uj'^ 



with a critical wavenumber kc = ^ 1 3l(] , which corresponds to a perturbation with 

small diffusive length scale. The phase diagram in the parameter space a-uu represented in 
Fig. [1] was numerically reproduced starting with an initial condition with kxw = (1, 1) 27r/3 
plus noise of amplitude 0.02. The system has been integrated using periodic boundary 
conditions in a square box of length 3. The mesh of Fig. [1] represents the final states for 
the corresponding values of a and cj, for i = 10 and e = 0.0083 (each square in the mesh is 
the result of an individual numerical integration). To the right of the parabola a = uu^ the 
traveling wave solution becomes unstable and gives rise to another structure with smaller 
wavelength associated with the diflFusive terms in Eq. ([5j). Some squares of the mesh still 
show the long wavelength solution to the right of the parabola, where it should be unstable. 
The reason is that, close to the stability limit, the unstable modes require a time greater 
than t = 10 to grow. 

The nonlasing solution, = 0, is linearly unstable for a > if < 0, and for a > —uu'^ 
if > (exhibiting again a large diffusive critical wavenumber kc = ^ 1) [31]. The 

numerical phase diagram of Fig. [2] confirms again the theoretical stability predictions and 
shows the appearance of a structure with wavenumber k ^ 1/v^ ^ 1 to the right of the 
stability limit given by the parabola a > —uJ^. The initial condition is Gaussian noise with 
amplitude 0.2, the final time is t = 10 and e = 0.0083. 
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In the next section we will analyze the dynamics in two representative points of the phase 
diagram: a = 0.75, uj = 0.5 (pattern with dispersive scale ^disp ^ 1), and a = 0.5, uj = 2 
(pattern with diffusive scale ^difr ^ V^)- 



III. NUMERICAL VALIDATION 



We numerically check the accuracy of the CSH equation ([5]) as a reduced dynamics of the 
Maxwell Bloch equations ([I])-([3]). The difference between both descriptions is computed as 



d = 





E{x,y,t) 




1 




P{x,y,t) 




1 




N{x,y,t) 








a 



(7) 



Symbols 1 1 • 1 1 denote the Euclidian norm on C^^ divided by ViV, where is the number of 
points of the discretized system. So, d is an average absolute error, that, according to the 
weakly nonlinear procedure applied to the MB equations to derive the CSH equation, has 
to behave as ~ 5^; see Eq. (|6j). 

There is a severe numerical difficulty in integrating Eq. ([5]) due to the presence of two 
spatial scales, 5disp ^ 1 and 5diff ^ V^, which are very different in the relevant limit £ ^ 
and should be simultaneously well resolved. We consider a one dimensional system, with 
periodic boundary conditions, in order to reduce the size of the computations and be able 
to use a greater system length than would be possible in higher dimensions. We let the 
system evolve until the difference d reaches a stationary value rf^. For the parameters 
used, a time t ^ 10 is enough to reach a stationary state, and the corresponding maximum 
integration time for the MB equations is t ^ 10/^^. Therefore, to check the asymptotic 
theoretical behaviour for 6: ^ 0, we need a large number of Fourier modes and long time. 
The CSH and MB equations are integrated in Fourier space using a fourth-order Runge- 
Kutta scheme, with 1024 Fourier modes, a = 6 = 1, time step dt = 0.01 {di = 2dte^)^ space 
step dx = 1/1024 [dx = rfx/(26:)], and s in the range between 0.0011 and 0.025. 

The initial condition for is filtered Gaussian noise of amplitude 1. Since Eq. ([5]) does 
not include spatial scales smaller than (5diff ~ ^/s^ the modes with wavenumber greater than 
k ^ 1/v^ are initially filtered. The initial condition for (£", P, A^) is obtained from the one 
for (j) using Eq. ([6]). In Fig. [3] we show the real and imaginary part of at t = and t = 5, 
for £ = 0.0011, a = 0.5 and uu = 2. 
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In order to calculate an average relative error, we divide by 5, the typical magnitude of 
the fields in the Maxwell Bloch equations (note that (j) is of order 1). In Fig. H] we plot d/e 
in log scale against t, for different values of s and for a = 0.5 and u = 2. 

We consider the stationary value of the relative error for long times {d/e)s^ and check 
the theoretical asymptotic behaviour {d/e)s ^ £• In Fig. [5] we plot {d/e)s against e for two 
points in parameter space: a = 0.75, uu = 0.5; and a = 0.5, uj = 2. In both cases, the linear 
behaviour {d/e)s ^ is confirmed. But the numerical result offers a new important figure: 
the slope. The slopes are 3.6 ± 0.1 for a = 0.75, uj = 0.5; and 18 ± 1 for a = 0.5, cj = 2. 
The increase of the slope in the second case is related to the fact that patterns with smaller 
length scales appear for a = 0.5, uu = 2 (see Fig. [1]). 

An experimental confirmation of the CSH equation would require to know a specific value 
of the appropriate distance to threshold for which the equation is valid. Let us suppose that 
the sought experimental confirmation has a maximum relative error of 10% and make the 
favorable assumptions that the Maxwell Bloch equations accurately describe the experiment 
and that the chosen parameters correspond to a simple pattern with characteristic length 
equal to (^disp ^ 1 as, for example, for a = 0.75, uj = 0.5. Then, using the slopes of Fig.[5l we 
can calculate that the distance to threshold should not exceed r — 1 = 0.003, and the pump 
must be tuned with a relative error smaller than 0.3%. For a = 0.5, = 2, where patterns 
with diffusive scale 5diff ^ arise, the situation is worse since the maximum distance to 
threshold is r — 1 = 0.0006, and the relative error of the pump should be smaller than 0.06%. 

IV. CONCLUSIONS 

We performed numerical integrations of the Maxwell Bloch equations and the correspond- 
ing CSH equation, for a class C laser. The CSH equation gives a simplified and reduced 
dynamics of the original Maxwell Bloch equations for small detuning. Comparing the results 
produced by both set of equations, we obtain an average relative error of the CSH equation 
solutions. The numerical results confirm the following theoretical prediction: {d/e)s ^ ^, 
where {d/e)s is the stationary relative error that is reached for long times (the small pa- 
rameter e is introduced in the deduction of the CSH equation and is directly related to 
the detuning and distance to threshold). Therefore, as expected, the CSH equation with 
terms of different asymptotic order [31] is the appropriate envelope equation to accurately 
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represent the behaviour of the Maxwell Bloch equations when s ^ 0. 

The numerical results also allow us to estimate the distance to threshold range for which 
the CSH equation holds. Assuming an average relative error of 10%, the maximum distance 
to threshold is between 0.003 and 0.0006 for the parameter values analyzed: a = 0.75, 
(jj = 0.5, and a = 0.5, = 2. The most restrictive value (0.0006) corresponds to the case 
when the resulting pattern has diffusive scales ((5diff ^ ^/s). Although it is not unfeasable 
to experimentally establish such small distance to threshold, it requires a fine tuning of 
the pump that is not usually available in standard lasers. Another important experimental 
difficulty is to obtain a wide enough beam for the patterns to develop. 

Finally, it is important to mention that, despite of the problems for setting up an accurate 
laser experiment in the CSH range, the numerical results on this paper confirming the validity 
of the CSH equation are interesting and valuable from the more general point of view of 
Pattern Formation. The CSH equation is an envelope equation and it is universal in the 
sense that its structure depends only on the kind of instability of the problem and not on 
the particular physical problem under consideration. 
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FIG. 1: Real part of (j), in a 2D system of size 3x3 with periodic boundary conditions, for different 
points in the region of parameter space a > and a; > 0. Each square represents an individual 
numerical integration of Eq. ^ with initial condition given by a traveling wave, with Ictw — 
(1, 1) 27r/3, plus noise of amplitude 0.02. The curve a = o;^ is the stability hmit of this kind of 
solution. The final time is f = 10 and e — 0.0083. The gray scale limit values are: black ~ —2.5 
and white ~ 2.5 
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FIG. 2: Real part of (j) for different points in the region of parameter space a < and > 0. 
Each square represents an individual numerical integration of Eq. ([5]) with initial condition given 
by the nonlasing solution, = 0, plus Gaussian noise of amplitude 0.2. The curve a — —uj^ is the 
stability limit of the zero solution. The final time is t = 10 and e — 0.0083. The gray scale values 
are: black — —2.5, gray — 0, and white — 2.5 
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FIG. 3: Real (continuous line) and imaginary (dashed line) parts of in a ID system of size 1 and 
periodic boundary conditions. Top: the initial condition given by filtered Gaussian noise. Bottom: 
final state for t = 5. Parameters are s = 0.0011, a = 0.5 and u = 2. 
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FIG. 4: Average relative error d/e in log scale against time i for different values of e. Prom top 
to bottom, £ = 0.025, 0.018, 0.013, 0.0088, 0.0063, 0.0044, 0.0031, 0.0022, 0.0016 and 0.0011. Param- 
eters are a = 0.5 and u = 2 
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FIG. 5: Stationary relative error {d/e)s against e. Plus symbols correspond to a = 0.75, u = 0.5 
(slope 3.6 ± 0.1); and asterisks to a = 0.5, to = 2 (slope 18 ± 1). 
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